SECANT ZETA FUNCTIONS 

MATILDE LALIN, FRANCIS RODRIGUE, AND MATHEW ROGERS 

Abstract. We study the series tpsiz) := X)^j sec{mrz)n^'' , and prove that it converges 
under mild restrictions on z and s. The function possesses a modular transformation 
property, which allows us to evaluate tps(z) explicitly at certain quadratic irrational 
^ I values of z. This supports our conjecture that T^~''tpk{y/j) £ Q whenever k and j are 

'~^- positive integers with k even. We conclude with some speculations on Bernoulli numbers. 
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1. Introduction 



Let C(s) denote the Riemann zeta function. It is well known that C(2n)7r~^'^ E Q for 
n > 1. Dirichlet L- functions and Clausen functions are modified versions of the Riemann 
1^ . zeta function, which also have nice properties at integer points [Lew91]. Berndt studied 

C^ I another interesting modification of the Riemann zeta function, namely the cotangent zeta 

function [Ber76]: 

CN ' He proved that (1.1) converges under mild restrictions on z and s, and he produced many 

S2^ ■ explicit formulas for Ck{z), when z is a quadratic irrational, and A; > 3 is an odd integer. 

^T) ', Consider the following examples: 

^' C3^^ =-— ^, C5(^/2) 



cn ■ \ 2 / 45^5' ' ' 945\/2' 



Berndt's work implies that \f] Cfc(-y/7)vr~ G Q whenever j is a positive integer that is not 
a perfect square, and fc > 3 is odd. A natural extension of that work, is to replace cot(2;) 
with one of the functions {tan(2;), csc(2;), sec(z)}. We can settle the tangent and cosecant 



C^ , cases via elementary trigonometric identities: 



oo 



Etanfvrnz) , , , , v-^ cscfvrnz) , , . 
V^ = c«(^) - 2cs(2z), Y. s = c-(^/2) - Cs(^), 



n° -^ — ' n 



but it is more challenging to understand the secant zeta function: 

, , , v^ secfTrnz 
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The main goal of this this paper is to prove formulas for specials values of "08 (-2) • In Section 
2 we prove that the sum converges absolutely if z is an irrational algebraic number and 
s > 2. In Section 4 we obtain results such as 

V^2(V2) = --, V2(V6) = — . 

These types of formulas exist because ipki^) obeys a modular transformation which we 
prove in Section 3 (see equation (3.8)). Furthermore, based on numerical experiments, we 
conjecture: 

Conjecture 1. Assume that k and j are positive integers, and that k is even. Then 
V'fc(Vj)vr^^ G Q. 

The results of Section 4 support this conjecture, even though there are still technical 
hurdles to constructing a complete proof. For instance, we prove that the conjecture holds 
for infinite subsequences of natural numbers. The rational numbers that appear are also 
interesting, and we speculate on their properties in the conclusion. 

2. Convergence 

Since sec(7rz) has poles at the half-integers, it follows that ips{z) is only well-defined if 
nz Z+ 2 for cLuy integer n. Thus, we exclude rational points with even denominators from 
the domain of ips{z)- If -z = p/q with q odd, then ips{p/q) reduces to linear combinations of 
Hurwitz zeta functions, and (1.2) converges for s > 1. Convergence questions become more 
complicated if z is irrational. Irrationality guarantees that | sec(7rnz)| 7^ cxo, but we still 
have to account for how often | sec(7rnz)| is large compared to n^. The Thue-Siegel-Roth 
Theorem gives that | sec(7rnz)| <C n^+^ when z is algebraic and irrational, and this proves 
that (1.2) converges for s > 2. The case when s = 2 requires a more subtle argument. 
We use a theorem of Worley to show that the set of n's where | sec(7rnz)| is large is sparse 
enough to ensure that (1.2) converges. We are grateful to Florian Luca for providing this 
part of the proof. In summary, we have the following theorem: 

Theorem 1. The series in (1.2) converges absolutely in the following cases: 

(1) When z = p/q with q odd and s > 1. 

(2) When z is algebraic irrational, and s > 2. 

(3) When z is algebraic irrational, and s = 2. 

Proof of Theorem 1, parts (1) and (2). Let 2; be a rational number with odd denominator 
in reduced form. It is easy to see that the set of real numbers {sec(n7rz)}neN is finite. Let 

M = max^gi^ I sec(n7rz)|. Then we have 

I sec(7rn2;)| M 



n' n° 

It follows easily from the Weierstrass M-test that (1.2) converges absolutely for s > 1. 
Now we prove the second part of the theorem. By elementary estimates 

|sec(7rn2;)| = |csc(7r(n2;- 1/2))| < | j r, (2.1) 

\i L/C T) '^n 

where kn is the integer which minimizes |nz — 2 — A;^! . Now appeal to the Thue-Siegel-Roth 
Theorem [Rot55]. In particular, for any algebraic irrational number a, and given e > 0, 
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there exists a constant C{a,e), such that 



a 
If we set a = z, then (2.1) becomes 

I sec(7rn2;)| ^ 

Therefore we have 



P 



> 



C{a,e) 



i2+e ■ 



(2.2) 



n \z 



2k„+l 
2n 



<Cn 



l+£ 



sec(n7r2;)| 1 
^ — - < 



s-l-e 



and this imphes that (1.2) converges absolutely for s > 2 + e. Since e is arbitrarily small 
the result follows. D 

In order to prove the third part of Theorem 1 , we require some background on continued 
fractions. Recall that any irrational number z can be represented as an infinite continued 
fraction 

z = [ao;ai,a2,...], 

and the convergents are given by 

Pi 



which satisfy 



[ao;ai,... ,a^] = — 



Pe. = aiPi_i+pi_2, 
qe = aeqe-i + qe-2. 



(2.3) 
(2.4) 



Convergents provide the best possible approximations to algebraic numbers among rational 
numbers with bounded denominators. In other words, if < q < q£, then 



P 

z 



> 



In addition 



1 



qiQe+i 



> 



Pe 



Pe 

z 

Qe 



> 



(.{qe+i +qe)' 



(2.5) 
(2.6) 



Now we state a weak version of a theorem due to Worley [WorSl, Thm. 1]. 

Theorem 2 (Worley). Let z he irrational, k> ^, andp/q be a rational approximation to 
z in reduced form for which 



P 



< 



i2' 



Then either p/q is a convergent pe/qe to z, or 

p _ api + bpi^i 
q aqg + 6g^_i ' 

where a and b are integers. 



\a\,\b\ < 2k, 



Now we can complete the proof of Theorem 1 . The following proof was kindly provided 
by Florian Luca. 
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Proof of Theorem 1, part (3). Let kn be the integer which minimizes \nz — 2 ~ kn\- Let 
Wz denote the set of integers where the quantity is large: 

(logn)^ 



Then 



W, 



n&Wz 



n G 



nz 



k-n 



> 



n 



' ^ ^' ^^n(log?i)2 



re' 



n=2 



which converges. 

Now assume that re Wz- Then 



1 + 2A;„ 



2re 



< 



(log re) 



re' 



Consider the convergents of z. Let i be such that g^_i < 2re < (7^ . By Theorem 2 there 
are at most O ((logg^)^) solutions to 



V_ 
2re 



< 



(logre)^ 



re' 



(2.7) 



with -p ^TL (i.e. consider all values of |a|, \h\ < 2k = 2(logg'£)^) 
From equations (2.5) and (2.6) we have 



p_ 
2re 



> 



Pe 



> 



1 



qe{qe+i + qe) 



Combining this with equation (2.4) implies 



nz 



l + 2kn 



r? > 



n" 



> 



qLi 



> 



{qe+i + qi) 8qe {qi{ai+i + 1) + qe-i) 8g|(a£+i + 2) 



Hence, if re W^, we find that 



sec(re7rz)| qj{ai+i + 2) ae+iq"^ 



<€. 



re' 



(il-i 



<C 



iti 



Combining the Thue-Siegel-Roth Theorem (equation (2.2)) with (2.6), implies that if z is 
algebraic 



1 



> 



qiqi+i 






> 



C{z,e) 



Thus qgj^i <^ Qi ■ This allows us to place an upper bound on a^+i: 

ai+i < < ql 



Putting everything together gives the bound 

I sec(re7rz)| 



< 



re' 



,l-4e' 



and as a result 



E|sec(re7rz)| v-^ (loga^)^ v^ 1 



n^VKz 



SECANT ZETA FUNCTIONS 5 

Since g^+i = a^+ig^ + qi^i > qi + Qe-i, we conclude that qi > Fg, where F^ denotes the 
^th Fibonacci number. Since the Fibonacci numbers grow exponentially, we have 

f^ , 1 + ^5 

qe > -^, where tp = — - — . 

Setting (f = (/?^~^ > 1, we finally obtain 

oo 



E\sec(mrz)\ v-^ 1 1 

Thus, it follows that (1.2) converges absolutely when s = 2. D 



3. A MODULAR TRANSFORMATION FOR ll^k{z) 
We begin by noting the following trivial properties of ipsiz)'- 

M-^) =Mz), (3.1) 

Mz + 2)=Mz), (3.2) 

2^->,(2z) =Vs(z) + Mz + 1). (3.3) 

The main goal of this section is to prove that ipkiz) also satisfies a modular transforma- 
tion. We start from the partial fractions decompositions of sec(7rx) and csc(7r3;), and then 
perform a convolution trick to obtain an expansion for sec(7rx) csc(7rx2;) (equation (3.5)). 
Differentiating with respect to x then leads to the transformation for ipk{z) (equation 
(3.8)). The third author used this method to rediscover the Newberger summation rule 
for Bessel functions [Rog05], [New82]: 

JMJ-.{y)-^,= E ^. (3.4) 

^ ' n=~oo 

Equation (3.4) follows from applying the convolution trick to partial fractions expansions 
for J^.(y) and J-xiv) [Lom68]. 

Lemma 1. Let X-4('^) denote the Legendre symbol modulo 4. Suppose that x and z are 
selected appropriately. ^ Then 

1 sr^ Y-A(n)csc(7mz/2) v-^ sec (vrnfl + 1/z)) 

7rcsc(7r2x)sec vrx =— + 8x> ^ ^ ^^ ) „ ' ' - 2zxy \ ^ 9 o • (3.5 

zx ^-^ n^ - 4x^ -^ n^ - z^x^ ^ ' 

n=\ n=l 

Proof. Recall the classical partial fractions expansions [GR94]: 

TT (seciTTx) - 1) =16x^ Yl n (n^-^Ix^) ' ^^'^^ 

n=l ' 

^csc(7rx)--=2xy \„ ' „ . (3.7) 

k=l 



We assume that z is irrational and algebraic, and x is selected so that the denominators in the sums 



are never zero. 
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Both sums converge uniformly. Multiplying the formulas together, expanding via partial 
fractions, and rearranging the order of summation, we have 

7r(sec(7rx) - 1) 7rcsc(7rzx) =32zx3 V — — ~ ^^ 

\ / n>l 

k>l 

=s,\^ y (-i)'^'x-4(n) / 1 i_ 



n{k'^ — z'^n'^/4:) \z'^n^ JAl — z'^x^ k? — z'^x'^ 

k>l 



Z— ,■ ri (n'2 — ZLt-S"! I Z— .■ i-2 _ 7 



vfc+1 



^ n (n^ — Ax^) \ ^-^ k"^ — z'^v? /A I 
^ A:2 - z2a;2 I Z^ ^ („2 _ 4^2/^2) I 

/C — _L N't — J. / 

By (3.6) and (3.7) this becomes 

vr (sec(7rx) — 1) ( vr csc(7rzx) I =32x'^ } -—r- ( vr esc {■Knz/2) ) 

y zx J '^ n^ (n^ — 4x^) \ nz J 

+ 2^^'^' E fc2(fe2_^2^2) (sec(^A:/z) - 1) . 

This reduces to (3.5) after several additional applications of (3.6) and (3.7). We can split 
up the sums because all of the individual terms converge absolutely. For instance, we can 
prove that 

^ (_l)fc+isec(7rfc/z) 

fc=i 

converges absolutely, by showing that the summand is <C \sec{'Kk/ z)\k~'^ , and then ap- 
plying Theorem 1. D 

Theorem 3. Let E^n. denote the Euler numbers, and Bm denote the Bernoulli numbers. 
Suppose that k € 2N. Then for appropriate choices of z: 

(1 + ^)'^ Vfe f T^ V (1 - ^)'" Vfc 



1 + zj \1- z 

r ■\k k ., X (3.8) 



k\ 

m=0 



^(2—1 - l)B^Ek-J^ [(1 + zr^-' - (1 - z)— 1] . 



Proof. Recall the Taylor series expansions of cosecant and secant: 

7rxcsc(7rx) = - > 2"' - 2)B„^- ^, sec vrx = V K ^ ^ 



m! ■' — ' m\ 

m=0 m=0 



Expand both sides of (3.5) in a Taylor series with respect to x. Comparing coefficients 
yields the following identity for k G 2N: 

M! J2 (2-- - DB^E,.^ ( ' ) .-- = -2^ f X-4(n)csc(W2) ^ ^._r^^ (, + ,/,) . 
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Finally, let z i— >• {1 + z) and z i— >• {1 — z), and subtract the two results to recover (3.8). The 
cosecant sums vanish because csc(7rn(l + z)/2) = csc(7rn(l — z)/2) whenever n is odd. D 

We conclude this subsection with a conjecture on unimodular polynomials. We call a 
polynomial unimodular if its zeros all lie on the unit circle. We have observed numerically 
that the polynomials in (3.8) have all of their zeros on the vertical line Jie{z) = 0. Since 
the linear fractional transformation z = {1 — x)/{l + x) maps the vertical line to the unit 
circle, we arrive at the following conjecture: 

Conjecture 2. We conjecture that the polynomial 

J2 2'"(2"^ - 2)B^Ek^m. Q) (x - xmi + xf-^, (3.9) 

m,=0 ^ ^ 

is unimodular when k G 2N. 

This new family of polynomials is closely related to the unimodular polynomials intro- 
duced in [LR12], [LS13], and [HZ13]. 



4. Special values of tpk{z) 

Throughout this section we assume that k is a positive even integer. Let SL2(Z) denote 
the set of 2 X 2 integer- valued matrices with determinant equal to one, and let PSL2(2;) = 
SL2(Z)/(I, — I). Recall that PSL2(Z) can be identified with the set of linear fractional 
transformations. The usual group action is 

a b \ az + b 



c d J cz + d' 

It is easy to see that multiplying matrices is equivalent to performing compositions of 
linear fractional transformations. Consider the following matrices in PSL2(Z): 

Equations (3.2) and (3.8) are equivalent to 

i:k{Az)=Mz), (4.1) 

V'fc(Sz)=(2z + l)^-Vfc(^) (4.2) 

+ ^ E (2""' - 1) B„^Ek-m Q) {z + 1)'-" [{2z + ir-^ - (2z + If-'^ 

m=0 ^ ^ 

Every matrix C G {A, B) has a factorization of the form C = A^^B^^A^^ . . . , so equations 
(4.1) and (4.2) together imply that there exists a z-linear relation between ipkiz), ipk{Cz), 
and TT^. 

Now we outline a strategy to obtain exact evaluations of ipk{z). First select a matrix 
C, and then find the linear relation between ipk{z), ipk{Cz), and vr . Next choose z so 
that ipkiz) = ipkiCz). For example, if z = 2j + ^/2j{2j + 1) in (4.2), then Bz = z - 4j 
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whenever j is a non-zero integer. With some work the equation collapses to 



^Pk[^/2j{2I 



+ r 






E(2 

m=0 



m— 1 



1) BmEk- 



1 + 



2j 
2j+l 



m—1 



2.7 
2j+l 



m—1 



1 + 



2i 

2i+l 



fc-1 



1 



2j 

2i+l 



fc-1 



(4.3) 



which holds for j G Z \ {0}. Similarly, if we take z = \2j + 1\ + y^2j{2j + 1), then we 
arrive at 



V^fc 1 + V2j(2j + 1) 



m=0 



/fc! 



^ (2--1 - 1) i?„i?,_ 



X 



1 + 



21 ) 



m—1 



22+1 
2i 



m—1 



1 + 



2i+l\ 

2j- ; 



fc-1 



1 



2j + l 
2i 



fc-1 



(4.4) 



which also holds for j G Z \ {0}. The right-hand sides of (4.3) and (4.4) are invariant 
under the Galois action -^/x i— )• — \/x, so both formulas are rational with respect to j. 
Specializing (4.3) at /c = 2 and k = A yields 



^2 (\/2j(2J + 1) 
i^A (\/2j(2j + 1) 



TT 



7b f + 46i + 6 
8.7 + 3 



■K 

180' 



These results support Conjecture 1. Further evidence for the conjecture is provided by 
combining (4.3), (4.4), and (3.3), to obtain identities like 

2 

i,2 (V8j(2j + 1)) = y, (4.5) 

whenever j £ Z. 

In general, it seems to be quite difficult to evaluate ipkiVJ) for arbitrary positive integers 
j. This is due to the fact that (4.1) and (4.2) restrict the available matrices to a subgroup 
of PSL2(Z). If there exists a matrix C G {A,B) which satisfies C^/J = y/j, then we can 
always evaluate tpk{y/j)- We can construct candidate matrices by solving Pell's equation: 



X2 - jy2 



If we choose C to be given by 



C 



X jY 
Y X 



then it follows that C^J = ^'^^x = ^' ^^^ det(C) = X^ - jY'^ = 1. Pell's equation 
has infinitely many integral solutions when j > 1 [HW38], so there are infinitely many 
choices of C. The main difficulty is to select appropriate values of X and Y so that C 
factors into products of A's and B's. It is not clear if such a selection is always possible. 
Notice that {A,B) C ro(2), so it follows PSL2(Z) ^ {A,B). 
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We conclude this subsection by noting that il^2iz) = for infinitely many irrational 
values of z. Setting n = — 3j in Proposition 1 below, yields 

= 0, (4.6) 

for any non-zero integer j. It is worth emphasizing that ipkiz) is highly discontinuous with 
respect to z, so these types of results are not surprising. 

Proposition 1. Suppose that j and n are integers, and n ^ 0. Then 




_2 



,U,m^L,^^, 



Proof. Setting A; = 2 in (4.2) yields 



Iterating (4.8) gives 



'ib2iBz) = ih2{z) + ^— -75 — - — . (4.8) 

^^ ' 2z + r^ ' {2z + lf 6 ^ ' 

, „ . 1 , . nz (3z2 + Anz + 2) vr^ , ^ 
Znz + 1 [2nz + \y b 

The derivation of (4.9) is best accomplished with the aid of a computer algebra system 

such as Mathematica, because significant telescoping occurs on the right. Now consider 
the matrix 



\^ 2n 4jn + 1 J 



and notice that Czq = zq, where 

zq 
Thus by (4.9) we have 

'4>2{zo) =V'2(Czo) 

=V2 (i?"(^0 + 2j)) 



2i(2in + 1) 



n 



1 n(zo + 2j) (3(^0 + 2i)2 + 4n(zo + 2j) + 2) ^2 



2n{zo + 2j) + 1 ^^^ ^^ (2n(zo + 2j) + 1)2 6 ' 

We complete the proof by solving for "02(^0) and simplifying. D 

5. Speculations and Conclusion 

Assume that k is a positive even integer. Euler gave the following expression for 
Bernoulli numbers: 

Bk = -2-^C{k). (5.1) 

Bernoulli numbers are interesting combinatorial objects, so it is natural to ask if the 
rational part of V'fc(Vj) a-lso has interesting properties. This is an obvious question because 
V'fc(O) = C{k). For instance, the von Staudt-Clausen theorem gives a complete description 
of the denominators of Bernoulli numbers: 

Bk= V - + Integer. (5.2) 

(p-i)lfc^ 
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Is there an analogue of (5.2) for the rational part of i^kiVJ)'^ As an example, consider 
V'fc(\/6) which can be calculated from equation (4.3). In order to eliminate some trivial 
factors, we define 

^' ■= 71 1^^'^^^^- 

By (4.3) we have 

The first few values are /32 = —8/3, /34 = 508/5, and /Sg = —64896/7, which all have 
denominators divisible only by primes where p — l\k (as hoped for). The first instance 
where this fails is /320- The denominator of /320 equals 5 • 7 • 11, but 7 — 1 = 6 does not 
divide 20. 
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